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SUMMARY 


Cosmos  1335  (1982-07A)  was  Launched  on  29  January  1982  into  a  near-circular 
orbit  of  inclination  74°,  and  decayed  on  5  April  1987.  The  orbit  has  been  deter¬ 
mined  from  observations  for  26  epochs  between  September  1985  and  May  1986,  a  time 
when  the  orbit  was  experiencing  the  effects  of  31:2  resonance  with  the  Earth's 
gravitational  field.  About  1400  observations  were  used,  the  most  numerous  being 
those  from  the  US  Navy  Navspasur  system,  and  the  most  accurate  those  from  the 
Hewitt  cameras  of  the  University  of  Aston  sited  at  Herstraonceux  and  at  Siding 
Spring  in  Australia.  The  average  orbital  accuracy  achieved  was  about  70  m  radial 
and  120  m  cross-track. 

Analysis  of  the  changes  in  inclination  and  eccentricity  at  resonance  has 
yielded  useful  values  for  six  lumped  harmonics  in  the  geopotential  of  order  31: 
the  two  lumped  harmonics  of  even  degree  had  an  average  accuracy  equivalent  to 
1.5  cm  in  geoid  height;  the  four  of  odd  degree  had  accuracies  between  2  and  7  cm. 
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1  INTRODUCTION 

The  satellite  1982-07A,  Cosmos  1335,  was  launched  on  29  January  1982  and 
decayed  on  5  April  1987.  It  is  believed1  to  have  been  an  octagonal  ellipsoid 
in  shape,  of  lengtli  about  1.8  m  and  diameter  about  1.5  m;  its  mass  was  about 
550  kg.  The  initial  perigee  and  apogee  heights  were  482  and  518  km  respectively, 
the  initial  orbital  period  was  94.63  minutes,  and  the  orbital  inclination  74.0°. 

In  November  1985  Cosmos  1335  passed  through  the  condition  of  31:2  reson¬ 
ance,  when  the  crack  over  the  Earth  repeats  every  two  days  after  31  revolutions. 
The  aim  of  this  Report  is  to  compute  accurate  orbits  from  observations  during 
the  time  when  the  31:2  resonance  with  the  geopotentlal  was  affecting  the  orbit, 
and  then  to  evaluate  lumped  geopotentlal  harmonics  of  order  31  from  the  changes 
they  produced  in  the  orbital  inclination  and  eccentricity.  The  orbit  was  deter¬ 
mined  between  September  1985  and  May  1986  from  radar  and  optical  observations 

2 

using  the  RAE  orbit  refinement  program  PROP,  in  the  PR0P6  version  . 

2  THE  OBSERVATIONS,  ORBITS  AND  OBSERVATIONAL  ACCURACY 

2 . 1  The  observations 

The  orbit  of  1982-07A  has  been  determined  at  26  epochs  from  1398  observa¬ 

tions  between  23  September  (985  and  7  May  ( 986.  Table  1  gives  the  number  of 
observations  accepted  for  each  orbit  and  their  sources. 

Observations  came  from  three  different  sources,  the  most  accurate  being 
those  from  the  University  of  Aston's  Hewitt  cameras  at  the  RGO  site  at 
Herstmonceux  (marked  H  in  Table  1)  and  at  the  Siding  Spring  site  in  Australia 
(marked  A).  These  observations  were  available  for  6  of  the  26  orbits  and 
usually  have  an  accuracy  of  a  few  seconds  of  arc. 

The  second  group  consisted  of  35  visual  observations  made  by  volunteer 
observers  who  report  to  the  University  of  Aston,  and  these  observations  usually 
have  accuracies  between  2  and  5  minutes  of  arc. 

The  third  and  largest  group  of  observations  comprised  Navspasur  observa¬ 
tions  kindly  supplied  by  the  US  Naval  Research  Laboratory,  usually  with 
accuracies  of  about  2  minutes  of  arc.  There  were  1317  of  these  observations 
available. 

2 . 2  The  orbits 

The  orbits  were  determined  at  approximately  8-day  intervals  from 
September  1985  to  May  1986  with  the  aid  of  the  RAE  orbit  refinement  program  in 
the  PR0P6  version,  and  the  orbital  elements  at  the  26  epochs  are  listed  in 


Table  1 


Sources  of  Che  observations  accepted  in  the  final  orbits 


Orbit 

number 

Source  of 

observation 

- 1 

1 

Total  J 

1 

Hewitt  camera 

H  '  A 

_ | _ 

Visual 

US  Navy 

1 

1 

l 

3 

61 

64  ! 

2 

i 

i 

2 

34 

36  | 

3 

i 

18 

18  I 

4 

36 

36  j 

5 

i 

55 

-  ! 

6 

46 

46  J 

7 

s  ; 

1 

42 

»  1 

8 

4  '  2 

6 

49 

61  i 

9 

i 

2 

46 

48  j 

10 

< 

48 

48 

11 

i 

37 

37 

12 

4 

56 

60 

13 

1 

43 

47 

14 

l 

74 

74 

15 

7  !  4 

7 

36 

54 

16 

l 

! 

53 

17 

I 

| 

50 

50 

18 

I 

80 

80 

19 

57 

57 

20 

;  4 

54 

58 

21 

: 

47 

47 

22 

i 

4 

62 

66 

23 

; 

60 

60 

24 

2 

50 

52 

2  5 

9  '  4 

4 

57 

74 

2b 

( 

66 

66 

- 

1 

Tota  1 

28  .  18 

35 

1317 

1398 
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Table  2,  with  the  standard  deviations  below  each  value.  The  epoch  for  each 
orbit  is  at  1)0  hours  on  the  day  indicated,  and  the  PROP  program  fits  the  mean 
anomaly  M  by  a  polynomial  of  the  form 

M  -  M0  +  :i,l  *  *  V  '  +  M*t"  +  M5t'’  ,  O) 

where  t  is  the  time  measured  from  epoch  and  the  number  of  M  coefficients 
used  depends  on  the  drag.  For  1982-07.\,  which  was  in  a  near-circular  orbit  at 
a  height  of  about  370  km,  were  the  only  coefficients  needed  on  8  orbits; 

14  needed  up  to  ,  and  4  required  M^-M^  . 

The  value  of  e  ,  the  parameter  indicating  the  measure  of  fit  of  the 

observations  to  the  orbit,  varied  between  0.43  and  0.83,  and  had  an  average 

value  of  0.60.  The  6  orbits  with  Hewitt  camera  observations  are  generally  more 

accurate  than  the  rest:  for  inclination  the  average  standard  deviation  for  the 
6  Hewitt  camera  orbits  is  0.0007°  and  tor  the  other  orbits  0.0011°;  for  the 
eccentricity  the  average  standard  deviations  are  0.00000b  and  Q.QQOOll,  respec¬ 
tively.  For  all  26  orbits  the  average  standard  deviation  in  inclination  is 
0.0010°,  which  is  equivalent  to  about  120  m  in  distance,  and  0. 000010  in  eccen¬ 
tricity,  equivalent  to  70  m  in  distance. 

Fig  l  shows  the  PROP  values  of  inclination  with  their  standard  deviations. 

The  values  of  eccentricity  are  plotted  in  polar  form  in  Fig  2.  This 
diagram  is  interesting  because  the  angular  motion  of  perigee  changes  from  circu¬ 
lation  through  360°  initiallv  to  libration  about  v  -  90°  later.  This  happens 
because  the  circular  locus  of  the  perigee  under  the  action  of  odd  zonal  harmonics 
is  converted  into  a  spiral  by  the  effects  of  air  drag  (and  possibly  resonance). 
The  locus  initially  passes  below  the  origin,  but  on  its  next  crossing  of  the 
e  sin  o  axis  it  is  above  the  origin. 

2-3  The  accuracy  of  the  observations 

The  residuals  of  the  observations  have  been  printed  out  with  the  ORES 
computer  program  ,  and  have  been  sent  to  the  observers.  The  rms  residuals  for 
observing  stations  witli  3  or  more  observations  accepted  in  the  orbit  determina¬ 
tion  are  given  in  Table  3.  Tne  PS  Navy  observations  from  station  29  are  geocen¬ 
tric,  and  if  they  were  given  in  the  same  form  as  the  lopocenlric  observations 
their  angular  rms  residuals  would  increase  by  a  factor  between  3  and  10.  In 
calculating  the  rms  lor  the  visual  observers,  residuals  greater  than  twice  the 
rms  have  been  omitted  (the  numbers  used  being  shown  in  brackets). 
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Table  3 


Residuals  for  observing  stations  with  3  or  more  observations 
accepted  in  the  final  orbit  determinations 


Number  of 

rms  residuals 

Station 

observat  ions 

Ra  nge 

Minutes  of  arc 

Accepted 

Re  jec ted 

km 

RA 

Total 

1 

US  Navy 

215 

9 

2.4 

2.6 

3.6 

2 

U  S  Na  vy 

156 

23 

2.8 

3.9 

4.8 

3 

US  Navy 

198 

18 

3.2 

3.2 

4.5 

5 

US  Navy 

162 

12 

3.8 

2.5 

4.5 

6 

US  Na  vy 

185 

12 

1.9 

2.6 

3.2 

29 

U  S  Na  vy 

401 

9 

0.5 

0.2* 

0.5* 

2265 

Fa  rnhtiiD 

5  (4) 

0 

3.4 

4.9 

6.0 

2414 

Bour.iemou  th 

3 

0 

4.2 

1  .7 

4.5 

2418 

Sunningdale 

7  (6) 

1 

2.7 

4.6 

5.3 

2420 

Willowbrae 

15  (12) 

2 

4.0 

3.4 

5.2 

2659 

Herstmonceux  3 
(Hewitt  camera) 

28 

0 

0.06 

0.07 

0.09 

9652 

Siding  Spring 
(Hewitt  camera) 

18 

0 

0.03 

0.04 

0.05 

*  geocentric 


The  residuals  in  Table  3  are  rather  larger  than  usual,  mainly  because  the 
orbit  was  much  closer  to  the  Earth  than  in  most  previous  orbit  determinations 
of  this  type,  the  height  (over  the  equator)  being  between  345  and  390  km.  Conse 
quently  the  satellite  could  often  only  be  seen  at  low  elevation,  where  observa¬ 
tions  tend  to  be  less  accurate;  this  especially  affects  the  US  Navy  observations 
When  the  satellite  was  observed  at  high  elevation,  it  crossed  the  sky  at  a  rapid 
angular  rate,  thus  aggravating  the  effects  of  timing  errors,  particularly  by 
visual  observers.  Also  the  satellite  '«s  small,  and  too  faint  to  be  ideal  for 
visual  and  photographic  observing. 

In  view  of  these  difficulties,  the  overall  rms  residuals  of  the  Hewitt 
cameras  are  remarkably  good,  being  5  seconds  of  arc  for  the  Herstmonceux  camera 
and  3  seconds  of  arc  for  the  Siding  Spring  camera.  Since  the  residuals  combine 
the  orbital  and  observational  errors,  the  observational  errors  of  the  Hewitt 
cameras  are  likely  to  be  smaller  than  their  rms  residuals,  and  2  seconds  of  aL„ 
would  be  an  accuracy  consistent  with  the  results. 


10 


3  THEORY  FOR  THU  RESONANCE  EFFECTS 

This  theory  has  often  been  given  in  detail  before  (tor  example  in  Ret  A), 
and  will  only  be  summarized  here.  The  longitude-dependent  part  ut  the  geopoten- 

c. 

tial  at  an  exterior  point  (r,:i,X)  is  written  as" 


V 

r 


1=2  m- 1 


(cos 


cos 


mA  +  SA  sin 
itm 


(2) 


where  r  is  the  distance  from  the  Earth's  centre,  6  is  co-latitude,  *  is 

longitude  (positive  to  the  east),  p  is  the  gravitational  -'onstam  for  the 

Earth  (398600  km  Vs  ),  R  is  the  Earth’s  equatori.il  radius  (h  178.1  km), 

P™  (cos  6)  is  the  associated  Legendre  function  of  i>-der  :a  and  degree  i  ,  and 

C  and  S„  are  the  normalized  tesseral  harmonic  c«»et t  ic lent s  of  which  only 
im  £m 

those  of  order  m  =  31  are  relevant  here.  The  normalizing  factor  N,  ^  is 
given  by 


N 


?,m 


2(21  +  l)(i  -  a)! 

( >  +  m  )  ! 


(  ’) 


The 
f  icients, 
of  order 


rate  of  change  in  inclination  i  caused 

C  and  SA  ,  near  R :a  resonance  ma 
£m  £  m 

)  as 


by  .i  relevant  pair  ut  cuet- 
be  written  (ignoring  term*. 


di 

dt 


n 


sin  i 


-  )  F„  0  (k  cos 
l mp  £  pq 


m) 


i  “m+  1 


■|jc 


(U 


where  F„  is  Allan's  normalized  inclination  function  ,  is  a  function 

£mp  £pq 

of  eccentricity  e  for  which  explicit  forms  have  boon  derived  by  Gooding, 
denotes  'real  part  of'  and  j  =  /-I  .  The  resonance  angle  t  is  del ined  by  the 
equa t ion 


t  -  i(  0  +  !!)  +  r  ( .  -  .  ) 


(V: 


where  is  the  argument  of  perigee,  M  the  rnea.i  anomaly,  tb.e  right  as« •eli¬ 

sion  of  Uie  node  and  v  the  sidereal  angle.  The  indices  v,  q,  k  and  p  in 
equation  (a)  are  integers,  with  y  taking  the  values  1,2,  1  ....  ind  q  the 
values  0,  *1,  ’2,  tlu  *  qua t  ions  linking  £  ,  m,  k  and  j>  are:  in  -  *  f  ; 

k  =  -  X  -  q  ;  2p  ■  ' 


He  re  B 


31  and  a 


and  we  shall 


l  l  rin.s 


-  1 1 

wni"-»  .ire  usually  dominant.  The  values  of  2  to  b».  ».k  .1  'u*-.c  he  such  that 

2,  >  m  and  (9.  -  k)  is  even.  The  successive  coefficients  which  arise  (for  given 

*■« 

,  and  q)  may  be  grouped  into  a  lumped  harmonic,  written  u~ 


2 


where  >.  increases  in  steps  of  2  from  its  minimum  perm: sarnie  value  2r  ,  and 

the  are  constant  coefficients,  with  Q  -  1  .  'H.e  vaiues  oi  the  «•  can 

c  2U  c 

be  obtained  from  equation  (4),  and  R.H.  Gooding  has  written  a  computer  program 
LKOr  foe  their  evaluation. 

It  has  often  been  shown  (for  example  in  Kef  7'  l»m'  .  In  the  ec«ia»l«  n  (4) 
for  di/dt  ,  only  the  q  =  0  terras  are  important  ii  e  1-  »m*il.  li  q  -  0 
and  =  1  ,  then  k  =  u  =  2  so  that  the  relevant  lumper*  .  i<  v  f*t  -i/d. 

are  Cjj*  and  •  Since  (2  -  k)  must  be  even,  tue  vi  uc?  of  >.  in  th. 

summations  in  equation  (h)  are  52,34,  .3$ ,  .  Inc  numerical  values  >»  e 

for  1982-G7A  are  given  in  section  4. 

The  rate  of  change  of  eccentricity  produced  by  a  r  •]  ,.r>  pair  of 

coefficients  C  and  S  near  B :a  resonance  me v  h-  w^-tten 
2  m  2  m 


de 
d  t 


!  R 
n  - 


F.  G,  , 

imp  x-pq/ 


)q  -  *5 (k  <-  q)e 


■K 


j^+1(C 


P.m 


jSt(r)  oxp  j(y$  -  qw) 


I  , 


(7) 


2 

where  terms  of  order  e^  have  again  been  ignored. 

For  near-c  i rcuia r  orbits,  the  only  important  terms  in  (7)  are  those  with 
q  =  1  and  q  =  -1  ,  because  all  the  others  are  multiplied  by  powers  of  e  - 
With  3=2  and  y  =  1  >  the  value  of  k(=ya-q)is  2  -  q  :  thus 
(q,k)  =  (1,1;  and  (-1,3)  are  the  important  terms.  since  -  *.)  aosi  be  even, 
the  values  of  2  arising  in  the  summation:*  of  equ>-ic."  <'  N  *  •*.»  •  -  .. 

Numerical  values  for  19b*. -u7A  >u«  given  in  s  cil^n  . 
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The  fitting 

Before  the  effect  r. 

ther  perln. 


.)  ,s. 


The 


1  i  '.a  t  i  on 


must  be  cleared  of 


•  )A  L  1  )I!S  . 
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and  lunisolar  gravitational  effects  have  been  removed  by  use  of  the  PROD  computer 
program8  with  integration  at  1-day  intervals.  The  perturbations  due  to  the 
harmonic  are  recorded  with  each  PROP  run  and  have  also  been  removed.  The  effect 
of  air  drag  was  removed  within  the  THROE  computer  program8,  assuming10  an  atmos¬ 
pheric  rotation  rate  A  of  1.1  rev/day.  Other  perturba t ions ,  such  as  solar 
radiation  pressure,  and  earth  and  ocean  tides,  have  been  ignored,  as  they  are 
expected  to  be  less  than  the  standard  deviations  of  the  values. 

Although  the  raw  values  of  inclination  (Fig  1)  look  unpromising  for 

analysis,  the  corrected  values,  shown  in  Fig  3,  are  much  more  regular,  and  a 

theoretical  curve  from  integration  of  equation  (4)  was  fitted  with  Lhe  aid  of 

9 

the  THROE  computer  program  .  The  values  of  the  lumped  harmonics  obtained  from 
the  fitting  are: 

9-0  2  9-0  2 

10  Cljp  =  1.4  ±  1.7  ,  10  s3i  =  11.8  ±  3.6  .  (8) 

The  measure  of  fit  c  had  a  value  of  0.95.  (e“  is  defined  as  the  sum  of 
squares  of  the  weighted  residuals,  divided  by  the  number  of  degrees  of  freedom.) 
None  of  the  standard  deviations  needed  to  be  relaxed,  but  the  last  three  values 
were  omitted  from  the  fitting,  because  they  were  a  long  way  past  the  resonance 
and  thus  contributed  little,  and  also  because  they  did  not  fit  well. 

Exact  resonance  ($  =  0)  occurred  on  18  November  1985,  and  Fig  4  shows  the 
variation  of  $  with  time,  which  is  close  to  a  straight  line  near  the  resonance. 

The  fitting  in  Fig  3  seems  entirely  satisfactory:  the  total  variation  at 
31:2  resonance  is,  as  usual,  quite  small,  but  the  orbits  are  accurate  enough  to 
define  the  variation  fairly  well. 

4 . 2  Discussion 

The  values  of  the  Q  coefficients  have  been  calculated  with  the  aid  of 
the  PROF  computer  program,  and  the  lumped  harmonics  may  be  expressed  in  terms  of 
the  individual  harmonic  coefficients  of  order  31  as  follows: 


C32,3l 
+■  0.086  C 


0.060 


42,31 


C34,31 
+  0.233 


-  0.503  C36)31  -  0.452  C.8>31 

C,,  +  0.242  C,,  ,,  +  0. 154 

44,31  46,31 


-  0.181  C40>31 

^48,31  + 


(9) 
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;  i 


The  equation  for  S^’  Is  similar,  with  C  replaced  by  S  throughout. 

[f  we  make  the  usual  assumption  that  the  C  are  of  order  10  ^ H"  , 

lm 

all  terms  In  equation  (9)  of  degree  >40  are  less  than  122  of  the  first  term. 
Ignoring  these  high-degree  terms,  we  can  rewrite  (9)  as 

C3°’2|l  +  0(0.15)  (  -  C32i2]  -  0.060  C34>31  -  0.503  C36/)1  -  0.452  C38  31  , 

.  (10) 

with  a  similar  equation  for  S  ,  where  the  term  0(0.15)  represents  a  conflation 
of  the  error  terms  and  is  small  beside  the  errors  in  equation  (8),  which  are 
greater  than  30%. 

There  have  hitherto  been  only  two  accurate  analyses  of  31:2  resonance' 1  W 
and  it  so  happens  that  the  Q  coefficients  for  one  of  these,  Samos  2  at  98° 
inclination12,  are  not  too  different  from  those  in  equation  (10).  For  Samos  2, 


0,2 

3 


'32,31 


0.257  C,.  ,,  -  0.186  C,,  -  0.355 

34 ,  3  1  36,3  1 


'38,31 


OD 


Thus  there  is  likely  to  be  some  similarity  between  the  numerical  values  of 

C?’i  and  1?!“  from  the  two  satellites,  unless  the  (34,31),  coefficients 

9-0  2 

happen  to  be  large.  The  values  from  Samos  2  were:  10  C, J  =  -2.9  1  1.2  ; 

9-0,2  3 

10  S3’  ■  9.0  t  2.2.  The  values  in  equation  (8)  are  as  close  as  would  be 

expected  to  those  from  Samos  2,  and  have  errors  about  50%  greater,  as  Is  also 

to  be  expected  because  the  drag  of  1982-07A  was  double  that  of  Samos  2. 

We  may  roughly  estimate  the  error  In  geold  height  Implied  by  the  standard 
deviations  a  In  equation  (8;  as  Ro/Q  ,  where  R  Is  the  Earth's  radius  and 

l  /  1,  \2/^ 

-  iT'l  2  /  2\  (  - 

Q  =  *0/*  )  j  •  Here  Q  =  1.144  ,  and  the  errors  in  geoid  height 

corresponding  to  the  a  in  equation  (8)  are  0.9  cm  and  2.1  cm  for  C  and  S 
respectively.  The  corresponding  accuracies  from  Samos  2  are  0.7  cm  and  1.3  cm. 

5  ANALYSIS  OF  ECCENTRICITY 

5 .  1  The  fitting 

Before  analysing  the  values  of  e  ,  we  need  to  remove  the  perturbations  due 
to  zonal  harmonics  and  lunisolar  gravitational  effects,  though  the  latter  prove 
to  be  very  small.  This  was  done  with  the  PROD  computer  program,  using  1-day 
integration  steps  (checked  by  re-running  with  0.5-day  steps). 
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The  removal  of  air  drag  effects  within  THROE  is  not  satisfactory  for 
1982-U7A,  because  the  air  density  model  within  THROE  is  spherically  symmetrical, 
whereas  the  real  atmosphere  has  a  ’daytime  bulge'  of  high  density.  In  the 
conditions  experienced  by  1982-U7A,  at  heights  near  350  km  with  low  solar 
activity,  the  maximum  daytime  density  is  4  times  greater  than  the  minimum  nLght- 
time  density,  according  to  the  COSPAR  International  Reference  Atmosphere  1972 
The  changes  in  eccentricity  produced  by  an  atmosphere  of  this  character  were 
calculated  by  the  method  given  in  the  Appendix  of  Ref  1L,  and  Fig  5  shows  the 
total  correction  to  e  ,  both  for  a  spin* rim  l  iy  symmetrical  atmosphere  and  for  an 
atmosphere  with  the  day-to-night  variation.  It  will  be  seen  the  two  are  com¬ 
pletely  different:  in  a  spherical  atmosphere,  e  is  always  reduced  by  drag,  so 
that  Ae  is  positive;  with  a  day-to-night  variation,  the  sign  of  Ae  depends  on 
the  position  of  perigee  (see  equation  (9)  of  Ref  14).  After  making  the  atmos¬ 
pheric  correction,  with  day-to-night  variation  included,  the  values  of  e  were 
fitted  using  THROE  with  drag  set  to  Zero  (that  is,  with  =  0). 

There  were  difficulties  in  fitting  the  points  after  MU)  46480:  but  this 
was  a  long  way  after  resonance,  as  Fig  4  shows,  and  it  was  decided  to  concentrate 
or.  the  17  central  values.  The  fitting,  with  (y,q)  -  (1,1)  and  (1,-1)  is  shown 
in  Fig  6,  and  the  values  of  lumped  harmonics  obtained  were: 


10 


9-  1 , 1 
31 


26  , 


9  -  - 1 ,3 
0  r  31 


-41  ♦.  10  , 


9- —1.3 

10  S 3J  =  -58  -  16 


(12) 


The  measure  of  fit  t  had  the  value  of  1.98  and  two  of  the  standard  deviations 
were  relaxed  by  a  factor  of  2,  as  shown  in  Fig  6. 

The  fitting  in  Fig  6  is  reasonably  satisfactory,  but  not  as  good  as  might 
have  been  hoped,  as  shown  by  the  rather  high  value  of  e  .  However,  values  of 
y.  between  1.5  and  2.0  are  quite  usual  in  fittings  of  e  :  with  Samos  2,  for 
example,  e  =  1.62  .  This  probably  happens  because  it  is  difficult  to  remove 
the  effects  of  drag  without  error. 

5 • 2  Discussion 

The  equations  for  the  lumped  harmonics,  with  the  Q  factors  calculated 
with  the  aid  of  PROF,  are 
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,  1 
'31 


C31  , 3 1  1  -922  C33,31  ”  0-757  C35,31  +  0,558  C37>31  +  1.052 


+  0.783  C41>31  +  0.169  C43>31  -  0.379  C45>31 


(13) 


r1-3 

°31 


C31 , 31  0-825  c33>3i  0.967  C35  3]  0.513  C3?>31  +  0.049 


+  0.445  C 


41 


,3?  +  0-571  c43,3i  +  °-/*52  C45>3]  +  ...  , 


.  (14) 

with  similar  equations  for  S  . 

The  error  in  geoid  height  implied  by  the  standard  deviations  in 
equations  (12)  may  be  assessed  by  the  method  of  section  4.2,  with  Q  =  2.283 
for  (q,k)  =  (1,1)  and  Q  =  1.563  for  (q,k)  =  (-1,3).  The  resulting  accuracies 
in  geoid  height  are  7  cm  for  C  and  2  cm  for  S  ,  with  (q,k)  =  (1,1);  and  4  cm 
for  C  and  7  cm  for  S  with  (q,k)  »  (-1,3).  The  average  is  5  cm,  and  this  is 
very  satisfactory  as  the  average  from  Samos  2,  with  drag  half  that  of  1982-07A, 
was  4  cm. 

6  CONCLUSIONS 

The  26  orbits  of  Cosmos  1335,  determined  from  US  Navy,  Hewitt  camera  and 
visual  observations  at  the  time  of  31:2  resonance,  have  proved  accurate  enough 
to  evaluate  even-degree  lumped  harmonics  of  order  31  with  an  average  accuracy 
equivalent  to  about  1.5  cm  in  geoid  height,  and  odd-degree  lumped  harmonics  with 
an  accuracy  between  2  and  7  cm.  These  results  are  for  a  new  inclination,  74°, 
and  should  make  a  valuable  contribution  to  the  determination  of  individual 
harmonic  coefficients  of  order  31. 
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